Ratchet-Like Solitonic Transport in Quantum Hall Bilayers 
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The pseudo-spin model for double layer quantum Hall system with total landau level filling factor 
V = 1 is discussed. Unlike the "traditional" one where interlayer voltage enters as static magnetic 
field along pseudo- spin hard axis, in our model we consider applied interlayer voltage as a frequency 
of precessing pseudo-magnetic field lying into the easy plane. It is shown that a Landau-Lifshitz 
equation for the considered pseudo magnetic system well describes existing experimental data. Be- 
sides that, the mentioned model predicts novel directed intra-layer transport phenomenon in the 
system: unidirectional intra-layer energy transport is realized due to interlayer voltage induced mo- 
tion of topological kinks. This effect could be observed experimentally detecting counter-propagating 
intra-layer inhomogeneous charge currents which are proportional to the interlayer voltage and total 
topological charge of the pseudo-spin system. 

PACS numbers: 73.43.Lp, 05.45.Yv 



In the recent years there exists a steady interest in 
quantum Hall bilayers with total Landau level filling fac- 
tor i> = I due to their anomalous transport and tun- 
neling properties 0. Quantum Hall bilayers consist of 
electrons confined in closely separated two dimensional 
semiconductor layers in high magnetic fields. In the ab- 
sence of interlayer voltage, each layer of the system has 
a filling factor 1/1=1/2 = 1/2. Because of the layers are 
identical, an electron in one layer could be identified with 
a fake spin up, while in the other layer with a pseudo- 
spin down. Therefore, the system can be phenomenolog- 
ically described via the pseudo-spin formalism . The z 
component of the overall pseudo-spin vector specifies the 
charge imbalance between the layers. It is clear that the 
system has a lower energy when the pseudo-spin points 
neither up nor down, but rather lies in the plane, re- 
flecting the fact that in the ground state electrons are 
equally distributed between the two layers. Therefore, 
within this formalism, double layer quantum Hall system 
is treated as an easy plane ferromagnet with a hard axis 
anisotropy and an electron tunneling between the layers 
corresponds to the spin flips in pseudo-spin language. 

This is a quite satisfactory model for the isolated dou- 
ble layer systems. However, problems start to arise as 
soon as one considers a real experimental situation with 
applied interlayer voltage and induced tunneling current. 
Traditionally, in the phenomenological Hamiltonian, in- 
terlayer dc voltage is interpreted as a constant magnetic 
field along z axis 0. Although this model correctly de- 
scribes the physics at low interlayer voltages it fails 
to describe experimentally observed current- voltage char- 
acteristics llj for large interlayer voltages. Besides that, 
in order to capture the essential physics of the system, 
one has to introduce effective damping mechanism in the 
formalism. 

Here we use a different interpretation of the interlayer 



voltage, comprehensive analysis of which can be found in 
Ref. Particularly, we treat the interlayer voltage as 
a circularly polarized oscillating magnetic field with the 
magnitude equal to the tunneling amplitude. In addi- 
tion, it is important that the driving magnetic field has 
the frequency proportional to the interlayer voltage it- 
self. In the absence of damping, the traditional model 
coincides to our's one in the rotating (with the frequency 
of oscillating field) reference frame. Nevertheless, after 
introducing the damping, the ground state of the system 
in these two models differ from each other (spins aligned 
in easy-plane in our model and tilted from easy-plane in 
the traditional one) and consequently the equations of 
motion are no more invariant with respect of changing 
the reference frame. 

In the equation of motion we choose Landau-Lifshitz 
damping term since it is most convenient for the sys- 
tems with easy-plane symmetries. The point is that it 
does not lead to the damping to a definite direction in 
easy-plane. Therefore, it is not needed to break the sym- 
metry "by hand" This would not be justified, since 
the only physically measurable quantity is the charge im- 
balance between the layers, i.e. z component of pseudo- 
spin, while all other components of pseudo spin are not 
physically measurable. As we shall demonstrate below, 
our model well describes experimental observations by 
Spielman et.al. 0. Moreover, this model leads to the 
prediction of directed solitonic transport phenomenon, 
which can be directly verified on the experiment. Note 
that, this novel effect appears as a natural consequence of 
the identification of interlayer voltage with an ac driving 
magnetic field. Indeed, as it had been suggested recently 
by Flach and co-workers in Ref. , ac driving force with 
a zero mean value may cause a directed transport in a 
nonlinear system. In addition, it was shown that the nec- 
essary condition for this phenomenon to appear is that 
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the well known Hamiltonian density: 




FIG. 1: Suggested experimental setup for observation of inho- 
mogeneous counter propagating intra-layer charge currents by 
applying an interlayer dc voltage on a quantum Hall bilayer. 



the system is characterized by nonzero total topological 
charge. As is well established, such topological picture 
may exist only if the ground state of a system is degener- 
ated, what is the case in quantum Hall bilayers according 
to our model. 

Here we present the symmetry analysis and numeri- 
cal simulations in order to show that the directed inhomo- 
geneous intra- layer current appears due to propagation of 
topological excitations in a quantum Hall bilayer. Real- 
istic experimental setup is suggested in order to observe 
this exotic phenomenon. 

The Model: The effective Hamiltonian density of our 
phenomenological model for double layer quantum Hall 
(pseudo) ferromagnet is given by: 



dx 



+ /3(r 



^SAS 



.}(1) 



where m(x,t) is a order parameter unit vector; m± = 
•nix ± ivriy {niz = v\ — describes local electrical charge 
imbalance between two layers); is the in-plane spin 
stiffness, /3 gives a hard axis anisotropy, Agyis is a tun- 
neling amplitude, w = eV/Ti includes interlayer voltage 
V , e is an electron charge and "c.c." indicates a complex 
conjugated term. Then Landau-Lifshitz equations of mo- 
tion, which conserve the length of a local spin density, can 
be presented as follows: 
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is the effective magnetic field and 7 is the Landau-Lifshitz 
dimensionless damping coefficient. 

It is easy to see that in the absence of damping our 
model coincides with the traditional one 0]. Indeed, by 
redefining the order parameter as to+ to+ exp(— jwf), 
the same equations of motion (|2Jl could be derived from 
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Therefore in case of zero damping the Hamiltonians 
and Q describe the same physics in different frames. 
However, for nonzero damping these two approaches are 
no more equivalent. 

Numerical Experiment and Analytical Solutions: To 
specify the problem, we consider a quasi-one dimen- 
sional quantum Hall double layer system with the pe- 
riodic boundary conditions (quantum Hall "ring" as is 
shown in Fig. 1). By applying the appropriate dis- 
cretization procedure drh{x,t) / dx — > TOj+i(t) — rfii{t) 
and d^rh{x,t)/ dx^ — > rni+i{t) + rhi-i{t) — 2rhi{t), the 
equations of motion Q reduce to a set of ?)N ordinary 
differential equations {N is a number of discretization 
points). We solve this set of equations using MATLAB. 
The results are presented in Fig. 2. One can see that the 
charge imbalance is very small even for very large uj-s 
(comparable with j3). In the case of small voltages, as is 
expected, z component of the order parameter is slightly 
tilted from the easy-plane (the upper panel in Fig. 2), 
while for the large voltages, which corresponds to the 
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FIG. 2: Numerical simulation results: pseudo-spin evolution 
for low (upper graph) and large (lower graph) interlayer volt- 
ages. Main plots and insets represent pseudo-spin x and z 
components, respectively. Note that rriz is small in both low 
and large voltage limits. In numerical experiments it is taken 
/3 = 10, AsAS = 0.002, 7 = 0.01, u = 0.1 (upper graph) and 
u) = Z (lower graph). 
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large driving frequencies, z component of the order pa- 
rameter oscillates around the zero value (the lower panel 
in Fig. 2). Such behavior for the large voltages follows 
from the fact that the driving frequency is proportional to 
the voltage. As is well known |6j a high frequency driving 
field, even it has a large amplitude, can cause only small 
oscillations in a system around the equilibrium point. 

These observations suggest to simplify further the 
equations of motion. Particularly, in the limit — > 0, 
one can rewrite Eq. (j^J as follows: 



dt 
dip 
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2AsAS sin(v3 + ivt) (5) 
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where the phase variable ip is defined from the relation 
TO_|_ = — exp(i</3). Reminding that describes a 
local charge imbalance between the layers, Eq. |SJ) could 
be interpreted as the charge continuity equation with a 
damping term, where 

dip 

Js = -2epE-^, Jtun = 2eA5^s sin((y3 + ujt) (7) 

are intra-layer current in each layer and interlayer tunnel- 
ing current, respectively. Further, substituting © into 
the continuity equation finally we obtain: 

^^Pl-^ + ^-^PEP^+^P^SAssHp+ut) = 0. (8) 

For the first time, without the gradient term, this equa- 
tion was suggested by Wen and Zee (see also Ref. Q). 
The heuristical considerations in Ref. [7| as well as our 
numerical simulations of Eq. ^ suggest that, there are 
completely different regimes of the pseudo-spin dynam- 
ics for the low {eV/fi < V8/3Asas) and large {eV/h > 
\/8f3AsAs) voltages. In particular, in the first case under 
the periodic boundary conditions the analytical solution 
of ||SJ) should be sought as homogeneous time rotations 
around z axis with the angular frequency lo = eV/h: 



p = (po- ^t-, 



(9) 



where 4>q is a constant quantity. Substituting lO into the 
motion equation one can define (f>o as follows: 



2AsAssin0o = ^7- 



(10) 



Moreover, from Q the expression for the tunneling cur- 
rent reads: 



Jtun = {e^imi- 



(11) 



For the case of large voltages the solution is sought 
as homogeneous small time oscillations around a definite 
direction in the easy plane as: 



p = A sin((/)o -I- ujt) , 



A< 1, 



(12) 



where the constants A and 0o must be defined pertur- 
batively by substituting (|12|l into the equation of motion 
© (see Refs. |1,0). As a result one obtains: 
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and the dc component of tunneling current is given by: 



Jtun = 2eAsAssin A sin 
eAsAsAsiiKpo = - 
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It is easy to see that Eqs. (|11|) and H14() qualitatively well 
describe the experimentally observed tunneling current 
voltage characteristics 0. Indeed, at the low voltages 
the current increases with voltage, while for the large 
voltages the tunneling current decreases as . 

In the consideration given above it was assumed that 
total topological charge of the system is zero. Conse- 
quently we obtain that the intra-layer current is zero for 
any voltages according to the definitions The situa- 
tion, however, drastically changes if nonzero topological 
charge is present in the system. 

Directed intra-layer transport: In order to apply the 
symmetry analysis we should write down the expression 
for the density of energy flux in the system (see e.g. Ref. 
l5|]): 



dip d<p 
ox at 



(15) 



and define topological charge of the system as follows: 



Q = p{+(X)) — p{—oo) /27r 



(16) 



In general, the topological charge Q may take any integer 
value. For simplicity in our numerical experiment we 
choose an initial topological charge as Q = 1. That is, 
for the discretized problem of N "spins" we take: 



= exp 2i7r(j - 1)/{N + 1) 



0, (17) 



and apply the periodic boundary conditions rn^ = m^"*"^ 
to the system. According to the general approach Q let 
us consider the following symmetry transformations: 



p - 



-p. 



(18) 



These symmetry transformations leave the topological 
charge H16() invariant but change the sign of the density 
of energy flux l|15|l . In this case the symmetry proper- 
ties of equations of motion becomes crucial. If they are 
invariant under the symmetry transformations H18|l then 
the averaged energy fiux in the system is zero, otherwise 
directed energy flux exists even in presence of noise only 
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FIG. 3: Distribution of the intra-layer local charge current 
versus time and distance in the case of nonzero initial topo- 
logical charge {Q = 1) and nonzero interlayer voltage. In the 
numerical experiments the periodic boundary conditions are 
used. 



[^. The equations of motion (|HJ| are not invariant with 
respect the symmetry transform (|18f) and therefore, the 
directed intra-layer energy transport is expected for any 
appUed nonzero interlayer voltage. Similarly, for zero 
voltage Eq. ^ is invariant under the symmetry trans- 
formation lpl?|l and thus averaged energy flux should be 
zero. 

The numerical analysis of both, the initial motion 
equation ^ and its reduced version |(SJ) completely agrees 
with the predictions of the symmetry analysis. Thus, in 
the presence of nonzero total topological charge in the 
system, directed energy transport should be observed. 
Moreover, the direction of the transport can be changed 
just by inverting the sign of the voltage. Besides that, 
we can further extend the analytical consideration noting 
that the equation of motion in the " moving" refer- 
ence frame tp ^ ip — ut is nothing but dc driven-damped 
sine-Gordon equation (see e.g. Refs. The role 

of the dc driven force plays the term / = AP'yuj. As is 
well known, sine-Gordon equation in the absence of driv- 
ing force supports solitary wave solutions with nonzero 
topological charge, these solutions often are termed as 
kinks and are given by the following expression: 



(f = A arctan <^ exp 



-^SAS 

PE 



(X - Xq) 



(19) 



Applied driving force results in the motion of that object, 
velocity of which will be proportional to the driving force 
(applied voltage). Thus, by increasing of the voltage, 
the intra-layer energy transport (and consequently the 
inhomogeneous charge current) will be increased. 

It should be especially mentioned that if one would use 
the Hamiltonian Q as a starting point of the analysis, 
one gets the result similar to Eq. |(SJ) but with uj = 0. 
Consequently, no directed motion is possible according 



to the traditional model. 

Conclusions: The driven-damped pseudo-magnetic 
model is elaborated in order to describe the dynamics 
in quantum Hall bilayers at the total filling factor v = 1. 
The symmetry analysis and numerical simulations are 
used in order to show that the directed transport exists in 
the system in the presence of nonzero topological charge. 
Realistic experimental setup is suggested for observing 
the suggested phenomenon. Initial nontrivial topological 
charge [like presented in Exp. (|17|l ] in the system could 
be realized in the laboratory experiments by application 
of a weak in-plane magnetic field along the double layer 
"ring" {x direction in Fig. 1), for which a commensu- 
rate pseudo-spin distribution appears. Then, switching 
off the in-plane field and applying an interlayer dc voltage 
it will be possible to observe the inhomogeneous counter 
propagating currents in each layer (see Fig. 3). It is 
obvious that the thermal fluctuations will decrease the 
topological charge in a double layer svstern (like it hap- 
pens in narrow superconducting channels [13|) and as a 
result the intra-layer current eventually should decrease 
as well. 



Acknowledgements: R. Kh. is indebted to Ramin 
Abolfath and Kieran Mullen for stimulating discussions 
and is supported by USA Civilian Research and Develop- 
ment Foundation award No GP2-2311-TB-02 and NATO 
reintegration grant No PST.EV.979337. L.T. acknowl- 
edges financial support from the Deutsche Forschungsge- 
meinschaft (DFG) under Bu 1107/2-2 (Emmy-Noether 
program) . 



[1] 
[2] 

[3] 

[4] 

[5] 

[6] 

[7] 
[8] 

[9] 

[10] 

[11] 
[12] 



khomeriki@hotmail.com 
lasha@tkm.physik.uni-karlsruhe.de 

I. B. Spielman et al., Phys. Rev. Lett., 84, 5808, (2000); 
K. Moon et al, Phys. Rev. B, 51, 5138, (1995); K. Yang 
et al, Phys. Rev. B, 54, 11644, (1996). 
Y. N. Joglekar and A. H. MacDonald, Phys. Rev. Lett., 
87, 196802, (2001). 

M. Abolfath, R. Khomeriki and K. Mullen, Phys. Rev. 
B, 69 (April 15 issue), (2004). 

S. Flach et al., Phys. Rev. Lett., 88, 184101, (2002); M. 
Salerno and Y. Zolotaryuk, Phys. Rev. E, 65, 056603, 
(2002). 

L. D. Landau and E. M. Lifshitz, Mechanics, Pergamon 
Press, Oxford (1976). 

X.-G. Wen, and A. Zee, Phys. Rev. B., 47, 2265 (1993). 
M. M. Fogler and F. Wilczek, Phys. Rev. Lett., 86, 1833, 
(2001). 

R. E. Eck, D. J. Scalapino and B. N. Taylor, Phys. Rev. 
Lett., 13, 15, (1964). 

B. A. Malomed, Phys. Rev. B., 43, 10197, (1991). 

O. M. Braun et al., Phys. Rev. B., 67, 066602, (2003). 

J. S. Langer and V. Ambegaokar, Phys. Rev. 164, 498, 

(1967) 



